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Abstract 

First-principles calculations within the Local Density Approximation (LDA) or Generalized Gra- 
dient Approximation (GGA), though very successful, are known to underestimate redox potentials, 
fH , such as those at which lithium intercalates in transition metal compounds. We argue that this inac- 
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O \ curacy is related to the lack of cancellation of electron self-interaction errors in LDA/GGA and can 

_i ■ be improved by using the DFT-I-C/ method with a self-consistent evaluation of the U parameter. 

^vq . We show that, using this approach, the experimental lithium intercalation voltages of a number of 
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(T^ ' transition metal compounds, including the olivine Lia;MP04 (M=Mn, Fe Co, Ni), layered Li2^M02 
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^^ ' (x =Co, Ni) and spinel-like Lia;M204 (M=Mn, Co), can be reproduced accurately. 
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I. INTRODUCTION 



Redox processes are relevant to many technological applications, including corrosion, 
fuel cells and rechargeable Li batteries, and the ability to study these processes from first 
principles is therefore crucial. The key to a redox reaction is the transfer of electrons 
from one species to another. When the redox electron is transferred between very distinct 
environments (e.g. metallic to ionic) the standard Local Density Approximation (LDA) and 
Generalized Gradient Approximation (GGA) lead to considerable errors in the calculated 
redox energies. We show in this paper that treating self-interaction with the DFT+U^^^ 
method gives considerably better agreement with experiment and thereby provides a tool to 
accurately predict redox potentials. 

In particular, we focus on the study of Li insertion in transition metal compounds using 
GGA and GGA+U. Transition metal (TM) compounds have attracted intense research as 
cathode materials for rechargeable Li batteries due to their ability to simultaneously absorb 
Li^ ions and electrons. In the discharge cycle of a rechargeable battery Li is oxidized on 
the anode side and inserted as Li+ + e~ in the TM compound that comprises the cathode. 
The energy of this reaction determines the oxidization/reduction potential at which the 
battery operates. It is the high redox potential of Li cells that makes them so desirable in 
applications where high energy density is required. 

First principles calculations have been used extensively to predict important proper- 
ties of Li-insertion materials such as the average potentia l^i^i^i'''i^i^i^°i^'^i^^i^^i^^ and potential 
profil o^^i^^ for Li insertion, phase stability^iiSiiS and Li diffusion^S. While this has led to 
considerable success in predicting the trends of Li insertion voltages^ and even new phasesi^, 
it has been noted that LDA or GGA can give relatively large errors for the average Li in- 
sertion potential^'^^. For example. Table |l] compares the experimental voltage for different 
structures with the one calculated in the GGA approximation and with computational de- 
tails discussed in section UTTl The Li insertion potential is consistently underpredicted by as 





LiNiOs/NiOa 


LiMn204/Mn204 


LiFeP04/FeP04 


GGA 


3.19 


3.18 


2.97 


exp. 


3.8522 


4.1523 


3.524 



TABLE I: Calculated and experimental redox couple voltage in Volt. 



much as 0.5 to l.OV. Similar results have been obtained with LDA^. 

Recently, we have shown that electron correlation plays an important role in predicting 
the phase diagram of the Li2.FeP04 system^^, for which LDA and GGA qualitatively fail. 
In this work we demonstrate that the DFT+f/ method also corrects the voltage error from 
LDA and GGA. In our approach, U is calculated self-consistently^^S,, thereby making this a 
"first-principles" approach to predict redox potentials with no adjustable parameters. 

We first present some background information on the specific Li insertions materials 
investigated and how the electrochemical reactions take place in a rechargeable lithium 
battery. We also discuss the details of the DFT+^ method and the self-consistent calculation 
of U. In Sec. IIIII we show the results of our approach, highlighting the improvement over 
GGA and the good agreement with experiment. 

II. MATERIALS AND METHODOLOGY 
A. Materials and Crystal structures 

As a representative set of Li-insertion compounds, we have selected several materials 
representing different environments for Li and TM ions, which are well characterized exper- 
imentally. 

The family LiMP04 of olivine structures (M=Mn, Fe, Co and Ni) is a promising candi- 
date for rechargeable Li-battery electrodes in large applications such as electric and hybrid 
vehicles^l. Olivine-type LiMP04 and the de-lithiated structure MPO4, have an orthorhom- 
bic unit cell with four Formula Units (FU) and space group Pnma (see Fig. Q). The olivine 
structure can be thought of as a distorted hexagonal close-packing of oxygen anions, with 
three types of cations occupying the interstitial sites: 1) Transition metals M in corner- 
sharing MOg octahedra which are nearly coplanar to form a distorted 2-d square lattice 
perpendicular to the a axis, 2) Lithium ions in edge-sharing LiOg octahedra aligned in par- 
allel chains along the b axis, and 3) P ions in tetrahedral PO4 groups connecting neighboring 
planes or arrays. It is believed that the PO4 groups hybridize less with the TM than an 
oxygen anion does in simple close-packed oxides, and hence leads to more localized 3d states 
on the TM than in an oxide. 

The layered LiM02 and spinel-like Li2^.M204 are more traditional cathode materials that 



have been thoroughly studied experimentally^^ and theoreticallj*^*^*i2i2S. They are both or- 
dered rock salts (see Figs. I21E1)- The layered structure can be envisioned as two interpen- 
etrating fee lattices, one consisting of oxygen, and the other consisting of alternating (111) 
planes of Li and TM. In the R3m space group the Li and the metal ions remain fixed in the 
ideal rock salt positions, but the whole (111) oxygen planes can relax in the [111] direction. 
The spinel-like structure Lia;M204 is so named because at x = 1 it has the same structure 
as the spinel mineral MgAl204. We shall refer to it as spinel even when x = 2. It can be 
envisioned as a fee oxygen sublattice, with TM in half of the octahedral oxygen interstices, 
and lithium either in part of the tetrahedral sites at x = 1 or in the octahedral sites not 
occupied by the TM ions at a; = 2^^. 

B. Relation between insertion voltage and total energies 

When Li is inserted into a TM-oxide, its charge is compensated by an electron absorbed 
from the external circuit. The insertion reaction is symbolized by the following equation: 

AxLi + Li^MOy ^=^ U^+^JAOy, (1) 

where MOy is the TM compound host material. Using thermodynamical arguments, it is 
possible to relate the voltage V of the cell to the lithium chemical potential (/xlj) on both 
sides of Eq. [T] in the cathode^S: 

.cathode .anode 

V{x) = - ^^-(-^ ^ ^'' . (2) 

F is the Faraday constant, and /i|^°'^^ is the chemical potential in the anode, or more gener- 
ally, the chemical potential of the Li source. 

The average voltage {V) for Li insertion between two composition limits, Lij-^MOy and 
Li^^jMOy, can be found by integrating Eqn. |2l (usually between x = and 1), and is de- 
termined by the free energy of the compounds at the composition limit^. Neglecting the 
entropic and PAV contributions^, {V) can simply be determined by computing the total 
energy of LiajjMOj,, Li^^^MOy and Li: 

,,., _ - [E{Li,,MOy) - E{Li,,MOy) - (x2 - xi)E{Li metal)] 

^^ ~ {x2-x,)F ■ ^^^ 

Typically Xi = and X2 = 1 are taken as composition limits, as in these cases no Li- vacancy 

disorder occurs. 



Experimentally, the voltage vs. lithium composition curve V{x) can be conveniently mea- 
sured for both the charging and the discharging processes. The corresponding curves differ in 
general because of the overcharge potential present in the circuit. We obtain the experimen- 
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numerically averaging the charge and discharge 



32I over the appropriate composition range. 



C. The DFT-I-U method 

The DFT+U method, developed in the IQQO'si^, is now a well-established model to 
deal with electron correlation in TM and rare earth compounds. The method combines the 
high efficiency of LDA/GGA, and an explicit treatment of correlation with a Hubbard like 
model for a subset of states in the system. To investigate whether the underestimation of 
the lithium intercalation voltage in LDA/GGA could be related to Coulombic on-site effects 
we carried out rotationally invariant DFT+f/^ calculations. The essence of the method can 
be summarized by the expression for the total energy 

Elda+u[p, n] = ^lda[p] + ^Hub[^] - Edc[n] = £'lda[p] + Eu[n] (4) 

where p denotes the charge density and n is the TM on-site 3d occupation matrix. For these 
states the Hubbard interaction term -Enub replaces the LDA energy contribution E^c- The U 
correction term Eu = -EHub — -S'dc is defined by Eq. HI Although E^^ is not uniquely defined, 
we have chosen the spherically averaged version^^i due to the considerations discussed in Ref. 



Edc{n) = Trn(Trr?, — 1) = — ^Trn(Trn — 1), (5) 

Eu{n) = ^^Tr (n(l - h)) = ^Tr (n(l - h)) , (6) 

where we have defined the effective interaction parameter Ucs = U — J, or simply U af- 
terwards. The calculated energies are insensitive to the J parameter at fixed f/gfr^ and we 
include it in Ues- 

D. Self-consistent Calculation of effective U 



We determine the U parameter using the method presented in Ref. |26| which we briefly 



outline below. This method is based on calculating the response in the occupation of TM 



states to a small perturbation of their local potential. 

We start from an LDA/GGA {U = 0) calculation as the reference point. Then a small 
perturbation 

2 
m=-2 

in the local d-orbital potential is exerted on metal site i, where P^ represents the projector 
on the d states manifold of ion i, and a is the amplitude of the potential shift applied to 
the d levels. This induces a change in the occupation number of ion i as well as other ions. 
Thus we can calculate directly the response matrices, 

^'' dai' ^''' dai' ^'^ 

which measure the variation of the d-manifold charge density ra;^, on ion j, produced by a 
potential shift at ion i. The subscript "0" denotes the bare response, calculated without 
self-consistency (the Kohn-Sham potential apart from dV is frozen at the value obtained in 
LDA/GGA before the perturbation), and corresponds to the response from an independant 
electron system, while Xji is the screened response (charge density and potential relaxed to 
reach self-consistency). The effective interaction parameter U is then obtained as 

U = [X,' - X-%- (8) 

This is a well-known result in linear response theory, in which the effective electron-electron 
interaction kernel is given as a difference among the interacting density response and the 
non-interacting one^. Since DFT is used, a finite contribution from the exchange-correlation 
potential is also included in the effective U. As we use the integrated quantity n^ to probe 
the responses, the calculated effective interaction is averaged over the ion in the same spirit 
as DFT+f/. The matrix in Eq. |H1 whose diagonal term defines the on-site Hubbard U, also 
contains non-diagonal terms corresponding to inter-site effective interactions in LDA/GGA. 
These are not used in the DFT+f/ model. This method to compute U is contains full account 
of the screening to the external perturbation operated by the electron-electron interactions. 
In fact the perburtation is applied in larger and larger supercells until convergence of calcu- 
lated U is reached. We also notice that the calculation of U is based on the use of the same 
occupancy matrices entering the DFT+f/ functional, guaranteeing full consistency with the 
calculation we perfornt^. 



III. DETAILS OF THE CALCULATION AND RESULTS 

Total energy calculations are performed for Ni, Mn, Co and Fe in the olivines, layered 
and spinel structures whenever experimental information on the voltage is available. For 
each system the total energy of the lithiated and delithiated state is calculated with GGA 
and GGA+f/, with the projector-augmented wave (PAW) method^^^ as implemented in 
the Vienna Ab-initio Simulation Package^I. The use of GGA over LDA has previously been 
shown to be essential for correctly reproducing magnetic interactions and possible Jahn- 
Teller distortions^. An energy cut-off of 500 eV and appropriate fc-point mesh were chosen 
so that the total ground state energy is converged to within 3meV per FU. All atoms and 
cell parameters of each structure are fully relaxed. Jahn- Teller distortions are allowed where 
the transition metal ions are Jahn- Teller active (Mn^"*" and Ni^"*" in our case) by explicitly 
breaking the symmetry of the unit cell. Our relaxed cells of layered LiNi02 and spinel 
Li2Mn204 agree well with the calculations in Ref. |39|on Jahn- Teller distorted systems using 
GGA. All calculations are performed with spin-polarization. As discussed later, the total 
energy of a given structure depends critically on the magnetic state of the metal ions, and 
high-spin states are favored by the DFT-|-^ scheme we used. The ordering of the spin 
on the ions in different magnetic structures (i.e. ferromagnetic, antiferromagnetic or more 
complicated ordering) results in difference in the total energy of the order 10-60 meV per 
formula unit. From the total energies, the average lithiation potential can be calculated 
through Eqn. El 

Table ITTl shows the self-consistently calculated effective U values for Mn, Fe, Co and Ni 
in different valence states and structures. For each structure, U is calculated for the low 
and high valence states respectively in a fully lithiated and de-lithiated structure. In all 
cases, except Ni^^/Ni^"*" in the layered structure, a higher valence state leads to a higher U. 
For the three cases (Mn^+/Co^'^/Ni'^+) for which we have a [/ in a close-packed (layered or 
spinel) oxides and in an olivine phosphate structure, U is higher for the olivine structure. 
This may be related to the fact that the TM-octahedra in the olivine are only corner sharing 
in two directions but separated from each other by phosphate groups in the third direction, 
leading to very narrow bandwidth and well localized TM-d states. For comparison we also 
list the U values calculated in Ref. !40! for TM monoxides MO (M = Mn, Fe, Co and Ni) in 
non spin-polarized state. Good agreement with LiMP04 is found except for Fe^"*". We note 



that in Ref 26 the U value of 4.3 eV for FeO was obtained with the same hnear response 

n 

approach, in good agreement with Ref I40I. So the difference between our results and those 



in Ref 



401 could be mainly due to different crystal environment. 





Mn2+ 


Mn3+ 


Mn4+ 


Fe2+ 


Fe3+ 


Co2+ 


Co3+ 


Co4+ 


Ni2+ 


Ni3+ 


Ni4+ 


Olivine 


3.92 


5.09 




3.71 


4.90 


5.05 


6.34 




5.26 


6.93 




Layered 














4.91 


5.37 




6.70 


6.04 


Spinel 




4.64 


5.04 








5.62 


6.17 








Monoxide^ 


3.6 






4.6 




5.0 






5.1 







TABLE II: Calculated U in eV. 

Figures E] and El show respectively the average Li insertion voltage as function of U in 
the olivine, and in the layered and spinel structure. The horizontal short line indicates the 
experimentally measured voltage. Three calculated points for each system are marked on 
the curve: the small open circles indicate respectively the voltage one would obtain using the 
calculated U for the most reduced and most oxidized TM-state in each structure (e.g. Fe^+ 
and Fe^"*" in LiFeP04). The large filled circle corresponds to the voltage for the averaged U. 
The results for each system are discussed in more detail below. 

A. Olivine structures Lia:MP04 (M=Mn, Fe, Co, Ni) 

According to neutron-diffraction experiments^*^ the magnetic ordering of LiMP04 is 
antiferromagnetic (AFM) within the approximately square lattice of metal ions for each of 
the above four TM. FeP04 is also found to have AFM magnetic ordering^. The results 
in Fig. m have been calculated with AFM spin configuration in both end members. The 
calculated and experimental cell parameters, as well as the electronic occupation of the TM 
ions are listed in Table UTTl 

Mn Both Mn^"*" and Mn^+ are high-spin ions in GGA and GGA+U calculations. At- 
tempts to constrain them to lower spin states lead to much higher energy. FM ordered 
magnetic structures are 10 - 30meV higher in energy than the AFM ordered mag- 
netic structure as U is varied. A strong collective Jahn- Teller distortion is observed in 
MnP04, where Mn^+ is in the high-spin t2„e^ state, in GGA(+[/). The experimental 
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a (A) 


b(A) 


c(A) 


V(A3) 


TM ion config. 


LiMnP04 


GGA 


10.55 


6.13 


4.78 


309.13 


+3 p2 


GGA+U 


10.62 


6.17 


4.80 


314.52 


+3 p2 


Exp.30 


10.44 


6.09 


4.75 


302.00 




MnP04 


GGA 


9.92 


6.01 


4.93 


293.92 


+3 pi 
hg^g 


GGA+U 


9.98 


6.07 


4.96 


300.47 


+3 pi 


Exp3 


9.69 


5.93 


4.78 


274.67 




LiFeP04 


GGA 


10.39 


6.04 


4.75 


298.09 


+4 2 
hgSg 


GGA+U 


10.42 


6.07 


4.76 


301.07 


+4 p2 
hgGg 


Exp^ 


10.33 


6.01 


4.69 


291.39 




FeP04 


GGA 


9.99 


5.93 


4.90 


290.28 


f3 p2 
hgGg 


GGA+U 


9.99 


5.88 


4.87 


286.07 


+3 p2 


Exp.-2I 


9.82 


5.79 


4.79 


272.36 




LiCoP04 


GGA 


10.30 


5.93 


4.75 


290.13 


+5 p2 
hg^g 


GGA+U 


10.33 


5.97 


4.76 


293.55 


+5 p2 
hg^g 


Exp.31 


10.20 


5.92 


4.70 


283.90 




C0PO4 


GGA 


9.71 


5.48 


4.59 


244.24 


hg 


GGA+U 


9.98 


5.78 


4.74 


273.42 


+4 p2 
hg^g 


Exp^ 


10.09 


5.85 


4.72 


278.66 




LiNiP04 


GGA 


10.09 


5.91 


4.74 


282.66 


+6 p2 
hg^g 


GGA+;7 


10.12 


5.90 


4.73 


282.42 


fd p2 
hg^g 


Exp>3 


10.03 


5.85 


4.68 


274.49 




NiP04 


GGA 


9.66 


5.72 


4.71 


260.25 


+6 pi 
hgSg 


GGA+U 


9.92 


5.82 


4.84 


279.43 


+6 pi 



TABLE III: Cell parameters of the olivine structures in the lithiated and de-lithiated states 
well as the corresponding electron configuration at the TM ions. 



, as 



voltage for the MnP04/LiMnP04 redox couple has been obtained from Ref. |30l The 
voltage predicted with GGA+U (4.04 V at f/ = {Uy[^-2+ + f/j^„3+)/2) is within a few 
percent of the experimental voltage (4.1 V), and in sharp contrast to the large error 
made by GGA (Vgga = 2.98 V). 



Fe Both Fe^^ and Fe^"*" are high-spin in GGA(+t/) calculations, and the AFM ordering is 
more stable than FM ordering. Using Up^2+ and Up^3+ we calculated a voltage of 3.39 
and 3.55 V respectively. The voltage calculated with the average U= 4.30 eV is 3.47 
V, which agrees very well with the experimentally measured value of 3.5 V^^. This is 
a substantial improvement over the GGA predicted value of 2.97 V. Previously, the 
localization of electrons induced by U was also shown to qualitatively affect the phase 
behavior in this systen>2^. 

Co In LiCoP04 Co^+ is stable in the high-spin t^gc'l state. In the delithiated C0PO4, Co^+ 
is stable as non spin-polarized with GGA, but more stable by several eV with GGA+f/ 
in the high spin t2„e^ configuration at the calculated U value of 6.34 eV. As shown 
in Table IIIII the cell parameters of C0PO4 calculated with non spin-polarized Co^^ 
in GGA is appreciably smaller than experimental values, while GGA usually slightly 
overestimates cell parameters. With GGA+U and high-spin Go^^ the calculated pa- 
rameters are close to experimental values. While there is only limited electrochemical 
data on this material^i, the predicted voltage of 4.73 V at f/average is within a few % 
of the result 4.8V established by Anime et. a^, compared to the poor GGA predic- 
tion of 3.70 V. The high voltage of this material makes it particularly attractive for 
high-energy density applications. 

Ni Though LiNiP04 has been synthesized, no Li can be removed from it 
electrochemically^. Hence the voltage is probably larger than 5V, the limit of most 
electrolyte systems. At a; = 1 Ni^+ is stable as high-spin tSe^. At x = Ni'^"'" occurs 
in the low spin state ^20^0 ^r both GGA and GGA+U , but the high spin state ^L^g i^ 
less unstable in GGA+U than in GGA. Note that low-spin Ni^+ is a weak Jahn- Teller 
ion, and no appreciable collective distortion is observed in our relaxed unit cell. With 
f/average, & voltage of 5.07 V is obtained, which is in agreement with the fact that no 
Li can be removed from this material. 

B. Layered U^M02 (M=Co, Ni) 

For the layered and spinel structures AFM spin ordering on transition metal ions is 
topologically frustrated, and their actual magnetic ground states are not always clear in 
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experiment. But as the energy associated with different magnetic orderings is small, the 
simple FM ordering is used in the following calculations. 

Co In LiCo02 Co'^'^ is stable in the non spin-polarized state for the calculated Uq^3+ = 
4.91eV. At X = 0, Co^~*" is almost degenerate in either non spin-polarized or spin- 
polarized ^20 in GGA, but more stable with spin-polarization in GGA+U at the cal- 
culated ^co4+ = 5.37eV. While GGA+U still improves the agreement of voltage with 
experiment^ over pure GGA, the error for this system is larger than in the other sys- 
tems we calculated. This might be related to the fact that the GGA result is already 
closer to experiment than for all other systems. 

Ni In LiNi02 Ni'^^ is most stable in the low-spin tggC^ state and is a weak Jahn- Teller 
ion. With GGA a distorted unit cell is found with the short and the long Ni-0 bond 
length being 1.92A and 2.13A, respectively, compared to experimental values of 1.91A 
and 2.14A2S, and a stabilization energy relative to an undistorted cell of only -2meV, 
within the range of numerical errors, compared to -11 meV in Ref.^. With GGA+U no 
appreciable distortion is observed. Experimentally there is no cooperative Jahn- Teller 
distortion in LiNi02 though the Ni-0 octahedra are locally Jahn- Teller distortecP^, 
suggesting a very small stabilization energy,consistent with both GGA and GGA+U 
results. At a; = 0, Ni^"^ is stable as a non spin-polarized ion. The GGA+U voltage 
value of 3.92V agrees well with the experimental average voltage of 3.85V—, and is 
substantially better than the GGA result of 3.19V. 

C. Spinel Li^M204 (M=Mn, Co) 

For the spinel Lia;Mn204 there are two distinct plateaus in the voltage profile, between 
< a; < 1 and 1 < x < 2, respectively. For < x < 1 Li enters tetrahedral sites, while 
the reaction from LiMn204 to Li2Mn204 occurs through a two-phase process whereby the 
LiMn204 phase with only tetrahedral Li disappears at the expense of the Li2Mn204 phase 
with all Li octahedral. Calculations were done for x = 0, 1 and 2 structures to get separate 
average voltage values for the two processes. For M = Co the < x < 1 reaction potential 
curve is difficult to obtain accurately in experiments. Therefore only the average voltage for 
the 1 < X < 2 reaction is shown in fig. El 
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Mn Both Mn^"*" and Mn^+ are high-spin. Mn^+ is a strong Jahn- Teller active ion. In GGA, 
the calculated Mn-0 short and long bond lengths 1.94A and 2.40A agree with Ref. 



391 : in GGA+U they become 1.96A and 2.32A, respectively. Experimental values are 
I.94A and 2.29A, respectively^^, showing that the good structural prediction of GGA 
is retained in GGA+U. Goexistence of distinct Mn''"'' and Mn^"*" is found in GGA+U 
in the LiM204 compound. The GGA+f/average results (4.19V and 2.97V respectively, 
for the first and second plateaus) is in excellent agreement with the experimentally 
measured values of 4.15V and 2.95V—. 

Co Like in the layered structure, Go^+ in Li2Go204 is non spin-polarized, and at x = 
Go^"*" is more stable as spin polarized tjg in GGA+U. The GGA+U voltage (3.56V at 
f/average= 4.84eV ) agrccs very well with experimental data available for the LiiGo204 
to Li2Go204 reaction (3.5V-3^). 

Note that in the x = 1 structure of the spinel materials Lij;M204 we find distinct M^"*" and 
M^+ ions in GGA+U instead of ions of intermediate valence. The same phenomenon was 
observed in the intermediate structures Lia,FeP04 of the iron phosphate^^. This is a direct 
consequence of the Eu correction term to the total energy in Eq. 0] which penalizes the 
non-integral occupation of the rf-orbitals. Such charge ordering is necessary for correctly 
predicting the < a; < 1 and 1 < a; < 2 average voltage values of Li^.Mn204 simultaneously, 
as well as the 1 < x < 2 voltage of Li2.Go204, and is not present in pure GGA unless 
localization is assisted by a strong polaronic contribution such as the Jahn- Teller distortion 
around Mn^+. 

IV. DISCUSSION 

Introduction of Goulombic on-site correlations in GGA through the GGA+U clearly 
improves predicted lithiation potentials considerably over the use of pure GGA (or LDA 
for that matter). The errors of GGA+U and pure GGA on all systems for which we have 
experimental data are summarized in Fig. IHl Pure GGA consistently underestimates the 
lithiation voltage, which is a measure of the energy lowering when Li is transferred from Li 
metal (the anodic reference) to a Li+ ion and electron in the TM oxide or phosphate. The 
contribution of the Li"*" ion to the reaction energy is largely electrostatic, and one would 



12 



expect this effect to be well captured in GGA or LDA. Hence, the large voltage error in 
LDA/GGA must arise from the electron transfer from Li metal to the TM cation. Since the 
voltage is always underestimated in LDA/GGA these approximations clearly penalize the 
energy of the electron on the TM, thereby lowering the reaction energy. It seems reasonable 
to attribute this to the poor treatment of electronic correlations in LDA/GGA. In metallic 
lithium the electron is affected by a small self-interaction in LDA/GGA as its charge density 
is delocalized. On the TM ion, however, the electron occupies a much more localized d-orbital 
and will experience a much larger self-interaction. The lack of cancellation between the self- 
interactions contributions to the energy, which are related to an improper description of 
the correlation effects in LDA/GGA, leads to a systematic error in the prediction of the 
redox potential. In the direction in which the electron is transferred from a delocalized to a 
localized state, the reaction energy is penalized (not negative enough), making the potential 
too small. The use of GGA+f/ allows for a better description of the electronic correlation 
and, by discouraging fractional occupations of the Kohn-Sham orbitals, removes the spurious 
self-interaction thus producing a much more accurate prediction of the redox voltage. While 
we demonstrate the GGA/LDA problem and improvement obtained with DFT+t/ on Li- 
insertion materials, we believe that a more accurate description of correlation effects within 
the DFT+U scheme is also necessary in the study of other redox processes in which electrons 
are transferred between states of different kind (e.g. catalysis of organic molecules on TM 
surfaces). In fact, as explained in Ref. |2y, a better description of the electronic correlation 
(which enforces the independence of the single electron energy eigenvalues of the partially 
occupied states on their occupation, thus leading to the elimination of the spurious self- 
interaction) is needed to reproduce the physical difference among the ionization potential 
and the electronic affinity (or the band gap in crystalline solids) which plays a very important 
role in the energetics of processes involving electron transfer. 

In our calculations high-spin TM ions are always energetically favored by GGA+f/ over 
low-spin or non spin-polarized states. In C0PO4 the non spin-polarized Co^"*" in GGA leads 
to cell parameters inconsistent with experiment. In GGA+f/ Go^+ becomes high spin, 
improving agreement with experiment. For the other systems the GGA and GGA+f/ cell 
parameters are rather close, though GGA+f/ seems to lead to volumes that are slightly too 
high. Jahn- Teller distortions predicted by GGA are also reproduced in GGA+f/ for Mn^+. 

In summary, we have shown that the under-estimation of the lithium intercalation volt- 
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age in LDA/GGA can be corrected by using GGA+f/ with a self-consistently calculated 
parameters U, without sacrificing properties that are already accurately predicted by GGA 
(e.g. Jahn- Teller effect, cell parameters, magnetic ordering). Voltages for most systems are 
predicted within a few % of experimental values. 

We believe that DFT+f/ will significantly improve the accuracy of voltage prediction for 
candidate materials can be predicted, and therefore enhance the capability of screening new 
materials for their ability to be good cathodes. 
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FIG. 1: The olivine structure with cation polyhedra. 
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FIG. 2: The layered structure with MOg octahedra and lithium atoms. 
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FIG. 3: The spinel-like structure when fully lithiated (x = 2, Li atoms taking octahedral positions) 
with MOe octahedra and lithium atoms. 
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FIG. 4: Voltage as a function of U for the LiMP04 materials in the olivine structure. The 
short horizontal lines on the curves indicate the experimental voltage of the each material (no 
experimental information is available for LiNiP04). The two small open circles on a curve represent 
the voltage for U calculated in the oxidized (delithiated) or reduced (lithiated) states. The big 
solid circle represents the voltage at the average of the two U values. 



18 



L-Co 




^ U(eV) "^ 



FIG. 5: Voltage as a function of U for the layered and spinel structures. Legend the same as in 
Fig. H 
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FIG. 6: Difference between calculated and experimental voltage [5-9], for GGA and GGA+U, at 
the calculated U of the oxidized (delithiated) and reduced (lithiated) states, respectively (l=layered, 
s=spinel). For the spinel structures two voltage values for the < a; < 1 and 1 < x < 2 plateaus 
are calculated separately. Olivine LiNiP04 is not shown here because the voltage is unknown. 
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